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Turbulent temperature fluctuations in the atmospheric 
boundary layer measured at 2m, 7m and 23m in Ris¢g, Denmark, 
were analyzed with particular emphasis placed on determining 
characteristics of the high frequency region of the spectra 
of these fluctuations. The shape of the high wave number 
one-dimensional temperature spectrum and an estimate of the 
Kolmogorov scalar constant were determined. Comparisons 
of high frequency spectral regions of temperature and 
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I. INTRODUCTION 


Several attempts over the past few years have been made 
to measure high frequency temperature variations in the 
atmospheric boundary layer. Sensors with a sufficiently 
high frequency response and spatial resolution have been 
devised to overcome the inherent problems previously 
encountered in trying to measure small scale, high frequency 
memperature tilwetuations. In addition to the sensitive 
sensors, low-noise electronic equipment has become available 
to record and analyze small scale fluctuations. 

The object of this research was to utilize an analog 
analysis scheme to determine the shape of the high wave 
number (small scale) one-dimensional temperature spectrum 


and to evaluate the scalar constant, K,', which appears in 


9 
Piesexpression for the -5/5 form of the temperature spectrum. 
In addition, simultaneous analyses of velocity fluctuations 
mere conducted in order to obtain the necessary information 
Miitrmately required to evaluate Ky'- 

A knowledge of the small scale temperature distribution 
tie turbulent field has a number of applications. It is 
necessary when considering the heat budget which is of 
utmost importance in the fields of oceanography and 
meteorology. Sensible heat flux can be determined by 
the covariance of measured temperature fluctuations and 


vertical components of the simultaneous wind velocity 


MuGtuatlons . 
10 





Temperature inhomogeneities in the atmosphere are 
one of the quantities associated with the scattering of 
acoustic waves and electromagnetic radiation. This 
scattering, which is a major concern in the fields of 
acoustics and optics, is often related to variations in 
the refractive index. 

Further applications of a knowledge of the small scale 
temperature distribution are found in the areas of physi- 


cal chemistry and thermodynamics. 
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ee reteOne LOA DEVELOPMENT 


eee VELOCITY SPECTRUM 

Small scale turbulence in a fluid is a feature of 
mecrest to the environmental scientist due to its effect 
on heat, momentum and energy transport. The one-dimensional 


velocity (kinetic energy) spectrum is given by 


u~ = f 6(k)dk (Gh) 
0 

where u? is the total variance of the turbulent velocity 

fluctuations and k is the radian wave number. The rate of 

energy dissipation in an isotropic, turbulent flow field is 


given by 
e=15v f k“#{k)dk (2) 
0 


where the function k-8(k), the dissipation spectrum, 
Statistically describes the rate of decay of turbulent 
kinetic energy to heat. v(em?.sec *) is the kinematic 
Viscosity of the fluid. 

Modern turbulence research had as its impetus the 
Soncept of a continuous cascade of kinetic energy from 
Honee scale anistropic velocity fluctuations to smaller 
scales where the motion became increasingly isotropic and 


at which the energy is dissipated as heat. Small scale 


PZ 





motions which are characterized by small time scales are 
Statistically independent of the relatively slow large 
scale turbulence derived from the mean flow. 

This concept was presented by the Russian mathematician, 
A. N. Kolmogorov in 1941. Kolmogorov proposed the idea of 
fecal 1So0tropy, 1.e€., at some point in the cascade of 
lotic cnergysmrom large Scales to small scales, the 
turbulent fluid motion becomes isotropic. This idea was 
described by Kolmogorov in two hypotheses: 

(1) If the Reynolds Number (ratio of inertial to 
viscous forces) is sufficiently large, then there exists 
a range of small scales in which statistically, the turbu- 
lence is basically steady-state and is dependent only on the 
rate of dissipation of kinetic energy, e[(cem/sec)”.sec “)], 
the wave number, K(cm *) and the kinematic viscosity, 


v(em*.sec!). Accordingly, dimensional analysis yields 


o(k) = (ev?)"/4 F(k/K,) (3) 
for the form of the one-dimensional energy spectrum and 
Z Seal 
K(k) = k 2(ev?) 4/4 (k/k,)? F(K/KQ) (4) 


for the form of the energy dissipation spectrum. Here, 


Ke = (e/y2y 1/4 is the Kolmogorov wave number, 1/k. 1s known 
as the Kolmogorov microscale. F (k/k.) 1S a universal func- 
Elomewhich as valid for all fields of turbulence provided 


the Reynolds Number is sufficiently large. 


a: 





(2) If the Reynolds Number is very large (on the order 
of nO and greater) then a range of scales exists wherein 
the energy cascades independent of viscosity. In this 
range of wave numbers, the decay of the flow energy is 
characterized by a -5/3 dependence upon wave number. This 
range in which the energy is transferred inertially is 
known as the inertial subrange and characterizes the turbu- 
itenece Scales too large to be influenced by viscosity but 
small enough to be separate from the mean flow. In the 


miertial Subrange, 
6(k) = Kte2/5 , 9/5 (5) 


where 4(k) re the one-dimensional energy spectrum 
[(cm/sec)*.sec] de SiSmAaeUnyVversal Constant. Previous 
memaearcners have reported values of K* ranging from 0.48 
to 0.69 although the lower values are considered to be 


more accurate. 


Poe TEMPERATURE SPECTRUM 

It 1s possible to apply the Kolmogorov ideas to 
turbulent fluctuations of passive scalars such as tempera- 
ture. Temperature is a passive scalar when it is trans- 
ported in a fluid but does not introduce buoyancy effects. 


If 6(C°) is the deviation of the temperature from its 
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mean, the one dimensional temperature fluctuation spectrum 


is defined as 


ep? = sf o, (k)dk (C°)* . (6) 
0 
ire rate of temperature dissipation in an isotropic flow 
mee ld 15S 


alk 


See a k°o, (k) dk (CON? Bee (7) 
0 


s) 

K is the molecular diffusivity and ko, (k) is the scalar 
dissipation spectrum which describes the distribution with 
wawe numbers the rate of decay of the quantity yO the 
memperature Variance. 

Again, according to Kolmogorov: 

(1) If the Reynolds number is sufficiently large there 
exists a range of small scale temperature fluctuations which 


meatistically are dependent only on €, v, €, (the tempera- 


s) 
ture dissipation rate) and k. Once again dimensional 


analysis yields 
by (k) = ege 9/4 v4 HCG,k/k,) (8) 


for the one-dimensional temperature spectrum and 


ko, (k) =" oh aati go 


2 
: (k/k.)" H(o,k/k.) (9) 


Is 





for the temperature dissipation spectrum. o is the 
Prandtl number (o = v/k), the ratio of kinematic viscosity 
to thermal diffusivity. Ina gas o is approximately l. 
For air at 20°C the value is near 0.7. 

(2) If the Reynolds Number is even larger, then the 
temperature spectrum is independent of wv and Kk and is 
determined only by e and Eg: Dimensional analysis 
fWGor7~sin, 1951] yields 


6. (k) = A ae eS (10) 


Ko! 1s supposedly a universal constant although previous 
researchers have reported values ranging from 0.2 to 2.0. 
In view of this range of reported values of Ky! ste 
was the intent of this research to determine the shape of 
one-dimensional temperature spectra at high wave numbers 
and to evaluate Ky! from the data analyzed. In this way, 


a larger data base could be established from which deter- 


minations for Kg! could be resolved. 
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IIIT. INSTRUMENTATION AND EXPERIMENTAL APPROACH 


A. BACKGROUND 

Sensors and electronic systems with the necessary 
sensitivity and response to measure high frequency, small 
nease atmospieric temperature filvetuations were developed 
in the early 1970's by National Electrolab Associates 
Mamated of Vancouver, British Columbia, in consultation 
with the Institute of Oceanography, University of British 
Columbia. Previous electronics allowed the temperature 
Spectrum to be measured to just beyond the inertial sub- 
range but not to sufficiently high wave numbers to 
Doedquivocally evaluate the area under the scalar dissipa- 


tion spectrum. 


B. TEMPERATURE MEASURING SYSTEM 
1. Sensor 

The sensor of temperature fluctuations utilized 
was a resistance thermometer consisting of a platinum wire 
0.25 um (0.00001 in.) in diameter and 0.30 mm in length. 
The wire has a time constant of less than 10 ps in a 4 m/sec 
air flow. With a detection current of 50 uamps, its 
Sensitivity is 0.063 mv/°C. The sensor is made from 
Wollaston wire consisting of a silver jacket about a 
Platinum Core. The outside diameter of the jacket is 


45 um (0.0018 in.). The jacket and core are prestressed 


Ny 





into a V-shape and then soldered onto the end of a standard 
hot wire anemometer probe (Figure 1). The sensor has a 
frequency response which is essentially flat to 5KHz. 

2. Bridge 

Utilize antarn platinum wire to Measure tempera cune 
fluctuations takes advantage of the fact that a slightly 
heated wire is basically a simple resistance thermometer. 

In order to sense a change in resistance in the wire due 
fomevariations in the ambient temperature, a small heating 
current is passed through the wire and the resulting voltage 
mecp Gan be measured. If the heating is kept low, the 
mecistance Changes due to cooling by wind are negligible 
compared to resistance changes caused by temperature fluc- 
meaerons, A bridge circuit as utilized to ensure optimum 
current flow through the platinum core. 

The bridge circuit is a modified Wheatstone bridge 
consisting of two 2.21K2 resistors, sensor, probe and cable, 
balance adjustment Peeiicc ore Aideecavaed Ol Sar litem: 1c cd 
and variable capacitors (Figure 2) paralleling the balance 
resistance are necessary to match the probe and cable 
@amacitance. The bridge amplifier is a differential 
amplifier with a fixed gain of 60dB and a frequency 


mesponse tlat to 1S50KHz. 
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fee KEE REMENTAL PROCEDURE 
ie LOea tle 

The data were collected at the meteorological 
field site operated by the Meteorological Group of the 
Danish Atomic Energy Commission, Research Establishment, 
Risé. The site features a 130M meteorological tower 
instrumented at seven levels. In addition, miscellaneous 
instrumentation can be mounted at arbitrary levels provided 
oy that the electronics is sufficient to drive the cables 
from the instrumentation to the junction boxes at the 
seven standard levels. The tower stands on the Ris¢ 
Peninsula 0.5 Km from Roskilde Fjord. Signals from tower 
Mmistrumentation are Carried by cable to a meteorological 
Station 50m downwind of the tower where both analog and 
digital recording equipment are housed. During the course 
of the experiments winds blew essentially from the west 
off the fjord. Data were not collected unless the wind 
was from off the water. 

The Ris# Peninsula is 6Km north of the ancient 
Danish town of Roskilde which in turn is 30Km approximately 
southeast of Copenhagen on the island of Sjaelland. 

Zee 2 Uap ment 
a. Meteorological Tower Arrangement 
Mean wind speed and direction, temperature 

and humidity were measured at 7, 23, 39, 56, 72 and 96m 


(Figure 3). Sonic anemometer measurements were made at 


Za 
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ween ana /Zm. —Fine Structure turbulence measurements 
were made at different heights during the course of the 
experiment. In general, Naval Postgraduate School (NPS) 
eaeapment was deployed at 2m, 7m or 23m. Oregon State 
University (OSU) equipment was always at 23m. University 
of California, San Diego (UCSD) equipment was placed very 
near the top of the mast at 123m. Danish Atomic Energy 
Commission (AEK) equipment was located at 2m. Each group 
deployed sensors for measuring turbulent temperature and 
velocity fluctuations. The main purpose of the experiment 
was to make simultaneous measurements of atmospheric 
temperature fluctuations by different systems, at the same 
level and at different levels. The temperature sensor in 
each case was a platinum wire. The diameter of the NPS 
were was 0.25 um, of OSU 0.625 um and of UCSD both 0.25 
and 0.625 um. By making these simultaneous measurements 
with similar sensors but different electronics, it was 
hoped (1) to find in what way they differed and (2) more 
mnpertantly, to resolve some of the discrepancies surround- 
mg the value of the scalar constant Ky'- 

The main comparisons made in this thesis are 
between the NPS and OSU systems. 

b. NPS Platinum Resistance Thermometer 

At the 2m height the sensor was mounted in an 

aluminum holder which could hold as many as four such 


PensOus. FOrethne data discussed in this thesis the 


ZS 





temperature sensor was placed 5mm below a Disa hot wire 
anemometer. The other two positions were left vacant. 

At 7m the same aluminum holder was used but a 
Thermosystems probe was used for the hot wire. This probe 
was vertically aligned and then rotated to be within Smm 
to the side of the temperature sensor. 

A somewhat similar arrangement was made at 
23m except it was more difficult to get the temperature 
probe close to the hot wire probe. There was a separation 
of approximately lcm. 

Gaon Tre Anemometer 

A variety of hot wire anemometers was uSed. 
They were Disa (AEK, UCSD), Thermosystems (NPS, OSU), 
and a new weet vane probe system (AEK) designed by Larsen 
and Busch (1974). 

Disa probes were mounted in the aluminum holder 
previously described whereas the other systems had indi- 
vidual clamps. | 

3. Recording Systems 
NPS, OSU and UCSD used analog tape recorders. 
AEK recorded digitally on magnetic tape. The NPS data 
analyzed for this thesis were recorded on an Ampex FR 1300 
tape recorder at 7 1/2 ips using FM electronics. At this 
emoced the 3dB point 15 at 2.5KHz. Each NPS temperature 
Signal was recorded both direct and differentiated (with 


maopeetato tame). Ihe object of the differentiation is 
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to increase the level of the signal at high frequencies 
above the noise fevel. The differentiator circuit had a 
gain of 6dB/octave over the differentiating band of 
mecaucnciese, Thise band could besvaried from DC to 500Hz, 
bOUOHz, 1500Hz, 2000Hz or 3000Hz. The frequency at which 
unity gain occurred could also be varied. Beyond the band 
Manitterentiatton, the circuit rolled off at 6dB/octave. 
No other filtering of the signals was done. NPS hot wire 
Signals were recorded directly without filtering or 
differentiation. 

The input level of the tape recorder was set to 


handle 1 volt rms signals. 


IS 






IV. ANALYSIS PROCEDURES 


A. SELECTION OF DATA 

Generally when concerned with the analysis of high 
frequency turbulent temperature and velocity signals, 
it is not necessary to examine long sections of data. 

It was determined, however, in the preliminary analysis 

of this data that in order to obtain characteristic 
spectra, it would be necessary to extend the record length 
to 10 minutes. This was due to varying levels of the 
Tecorded signals, necessitating the utilization of a 
broadband amplifier to maintain the signal level at an 
almost constant level. 

The primary considerations for selecting a section of 
data for analysis were first that the temperature and 
velocity fluctuations be considered typical of the record 
and secondly, if possible, correspond to sections that were 
teins analyzed at Oregon State University for future conm- 
parison of results (Figure 4). By trying to adhere to 
these considerations, the following data sections were 
analyzed: a10minute length of record from Run 5 commenc- 
moe at 1537 local Ris® time on 30 August 1974, and a 10 
minute length of record from Run 6 commencing at 1657 on 


30 August 1974 (Table I). 
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B. ANALOG EQUIPMENT AND PROCEDURES 

The analog analysis scheme consisted of utilizing the 
following state-of-the-art equipment to obtain spectra of 
Emewselected sections (Figure 5): 

A Honeywell Model Ninety-Six Magnetic Tape Recorder/ 
Reproducer System was used to playback the recorded 
Signals. This device is capable of recording/reproducing 
any combination of multi-channel direct (analog) and FM 
data at nine selectable servo-controlled tape speeds 
maem lo/ homunrough 240 amches per second (ips). Its 
Standard magnetic assemblies are IRIG compatible. It is 
characterized by a tape speed accuracy of 0.1% and a 
bandwidth of 0-2.5KHz (within 1dB) at a signal-to-noise 
ratio of 50dB when operated at 7-1/2 ips. 

The playback signais from the reproducing system were 
fed to a Preston 8300 XWB Amplifier with a bandwidth of 
mowKHiz “and Selectable gains of 1, 2, 5, 10, 20 and 50. 

The output from the amplifier was fed into the input 
of a Federal Scientific UA-500 Ubiquitous Spectrum Analyzer. 
The selected analysis range was 0-2KHz at a sample rate 
of 6000hz. The analyzer has a 3dB bandwidth of 6hz and 
eperequency accuracy of +0.2% of the full analysis range. 

The input signal to the analyzer was monitored by a 
Hewlett-Packard 3400A RMS Responding Voltmeter and a 


Tecktronix Type 502A Dual-Beam Oscilloscope. 
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The spectral output was visually monitored on a 
HP 1300A X-Y Display and plotted on a Mosely 7035B X-Y 


Solid State Recorder produced by Hewlett-Packard. 


C. DATA INTERPRETATION 

Theoretical work almost always expresses turbulence 
results in terms of wave number while actual turbulence 
Measurements are spectrally computed as functions of 
frequency. These two approaches are related by Taylor's 


hypothesis 
k = —— Gil) 


where k is radian wave number cian f is frequency (Hz) 
and U is mean wind speed (cm/sec). Equation 11 implies 
that turbulence with scales near 1/k are "frozen" during 
time period these scales are swept past the sensor. The 
maximum frequency analyzed was 2KHz at a mean wind speed 
of S m/sec which indicated that a maximum wave number of 
25cm 1 VOMmdmiave tO mMeanhesorved , | [his requires that the 
sensor does not exceed 0.4mm in length and have a fre- 
quency response essentially flat to 2KHz. Both of these 
requirements were met. 

Since Eq. 11 is linear, the shape of a wave number 
spectrum is the same as that of a frequency spectrum. 


mrcduency Spectra and wave number spectra are related by 


Suk 





EO (1) Saeko Ch). (12) 


The inertial subrange expression (eq. 5) in terms of f 
becomes 


27 -2/3 


fey (ey reels pees 
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The rate of kinetic energy dissipation expression (eq. 2) 


becomes 
ee, 
—é= 15v im PsP an (56 cle (14) 
0 


The resulting expressions for the.temperature spectrum in 


terms of f are 


,(f) = (oe Toa eo (15) 
and 
oi 2 ery 
eae 6K (F- > is >, (f) df . (16) 


Equations 13, 14, 15 and 16 were the expressions used 


gd en Oa 


to evaluate ce, Eg 9 
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PeeeevELOCITY SPECTRUM RESULTS 
meee Calculatiomof ¢€ 
a. Indirect Method 
An indirect estimate of € was made from the 
friction velocity and the logarithmic wind profile 


eamecure 6). Thus, 


24%, 
Zo 


U, = KU/2N (iz) 
where K = 0.4 is the von Karman constant and Z) the reugh- 


Besos length. The relation 
ne 


was then used to evaluate e. 
b. Inertial Range Method 
Once the spectrum $¢(f) was determined, € was 
estimated by choosing values of $¢(f) from the best straight 
imimerin the frequency region in which a -5/3 form was 
Sesenved and then substituting into Eq. 13. A value of 
0.5 was assumed for the constant K! based on results 
presented by Boston (1970) and Williams (1974). 
Gee Dance ase tilod 
The kinetic energy dissipation rate was 


estimated directly by summing of the frequency range 
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Peon lOHz to 2hhiz the scales contributing to viscous dissi- 
pation according to Eq. 14. 

table Il summarizes the results obtained 
from these three methods. All values of e€ determined by 
the three methods were felt to be representative of the 
actual value. Therefore an average value for each run 
at a given height was utilized in the further analysis 
of the temperature spectrum. 

Z. Normalization of Spectral Results 
a. Velloemty Spectra 

The velocity spectra were normalized according 
to Eq. 3 so that a comparison could be made between dif- 
Zeremit runs and Previous experiments. Figures 7 and 8 
aacenoulalized velocity spectra from Runs 5 and 6 at the 
2m height. 

These figures show that the spectra nearly 
overlap with a characteristic break from the -5/3 region 
near log k/k. = -]., 

b. Energy Dissipation Spectra 

The energy dissipation spectra were normalized 
according to Eq. 4. A linear plot was used to display 
results (Figure 9). Maximum dissipation occurred near 
eae Value of k/k. = 0.1 which is in agreement with previous 


mesearch [Boston, 1970]. 
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Figure 7. 
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Bo TEMPERATURE SPECTRUM RESULTS 


ieacalrculatiaons Of Eg 


€, Values were estimated from Eq. 16 and the 


8 
mesultes are included in Table III. 


Ze Caveulations of Ka! 


The values of Kg! were determined from Eq. 15. 
Values of o,(f) were chosen from the best straight line 
in the frequency region in which a -5/3 form was observed. 
The value of e€ used was the average value for a given run 
fmanoarticwlar height. Since € enters into Eq. 15 to the 
-1/3 power, errors in estimating e€ do not seriously affect 
Ky'. 

The results of the Kg! calculations are shown in 
Table III. An average value of 0.88 was obtained which is 
in quite good agreement with 0.81 by Boston (1970). It 
Was noted that although no inherent differences in tempera- 
ture spectral shape were observed between OSU and NPS 
sensor results, the values of Ky! determined at the 23m 
meweht from the OSU sensors were very near the values 
reported by Williams (1974) at OSU. 

Se Normlalizatpton Of Spectral Results 
a. Temperature Spectra 
The temperature spectra were normalized 

meconding to Eq. 8. The temperature spectra from each 


run at each height were normalized and plotted (Figures 


10-15). A composite spectrum was obtained by superimposing 
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fisenie wspectra on a single plot (Figure 16). The break 
from the -5/3 region occurs just beyond log k/k. = -l. 
ners. spectra Show the shape of the one-dimensional 
temperature spectra in air beyond the -5/3 region which 
was one of the primary objectives of this thesis. The 
slope of the spectrum beyond the -5/3 break tended to be 
-32/5 which agrees well with Heisenberg (1948) who demon- 
strated that a slope of -7 should exist for very large 
values of k. With the exception of the temperature 
eee@era at 2m, there 1s very little scatter in the data 
at high wave numbers. 

PEGutcon!/ 1S aecompositre plot of Characters tre 
temperature and velocity spectra. It is noted that there 
1s very little difference between the two spectra but the 
temperature spectrum extends to slightly higher wave 
numbers than the velocity spectrum. A -32/5 slope is 
Observed for both temperature and velocity. 

Dem lempera tune Dissipation svectra 

The temperature dissipation spectra were 
fiomiaiized according to Eq. 9 and again plotted linearly. 
Figure 18 shows a typical normalized temperature dissipa- 
tion spectrum. There is very little difference between 
the temperature and velocity dissipation spectra. The 
normalized shapes are almost identical with the temperature 
Spectrum dropping off slightly faster than the velocity 
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Figure 16. Composite Normalized Temperature Spectra 
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VI. SUMMARY AND CONCLUSIONS 


The two primary objectives of this research were to 
determine the shape of the high wave number one-dimensional 
temperature spectrum and the evaluation of the Kolmogorov 
Scalar constant Ky! 

As seen in the results section, this investigator has 
determined that the high wave number one-dimensional 
temperature spectrum falls off with a slope of -32/5 beyond 
mie e=5/5 break Which occurs slightly beyond log k/k. = Jae 
Also the determination of the Kolmogorov scalar constant, 
Kay's woe bDemOeSeuls CoOMlsadered sto be Slenificant not omy 
in absolute value but also that it increases the data 
base of high frequency temperature fluctuation analysis. 
Further, the criticism that previous evaluations were based 
on records that were too short is removed. It was also 
observed (Table III) that. K,! increased in value both with 
height and Reynolds Number. 

As an additional result of this research, it is con- 
Cluded that the spectra of temperature and velocity at 
high wave numbers are very Similar with the temperature 
Spectra extending to slightly higher wave numbers than 
micmVvelOcilty Spectra but that they both fall off at the 


Same rate. 
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